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A simple method to deal with four dimensional Hamilton-Jacobi equation for null 
hypersurfaces is introduced. This method allows to find simple geometrical condi- 
tions which give rise to the failure of the WKB approximation on curved spacetimes. 
y^Q I The relation between such failure, extreme blackholes and the Cosmic Censor hy- 

pothesis is briefly discussed. 
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1 Introduction 



Hamilton-Jacobi (henceforth, HJ) equation is one of the most deeply analyzed 
equation of theoretical physics because of its application in many fields of great 
interest such as theoretical optics, quantum mechanics, Hamiltonian dynamics. 
In gravitational physics many important applications are based on H J equation 
as, for instance, gravitational lensing. HJ equation also encodes the light- 
cones structure of space-time and its analysis is also fruitful to disclose the 
asymptotic symmetries of space-time itself. In blackholes physics also, since 
the Carter's works [6] [7], HJ equation played a prominent role in clarifying 
many important local as well as global properties of such interesting objects. A 
seemingly new method to deal with the curved four dimensional HJ equation 
seems to be able to disclose interesting relations between complex structures, 
BPS states in gravity and the breaking of the WKB approximation on curved 
space-time. In particular, this method discloses in a clear and simple way the 
relation between the impossibility to adopt the geometric optics approximation 
on certain curved backgrounds and the presence of naked singularities. 
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2 The method 



In the limit of small wave-lenghts, beams of null geodesies describing, for 
example, electromagnetic or gravitational waves can be analyzed through the 
WKB approximation. The rapidly oscillating factor multiplying the slowly 
varying amplitude of the wave is usually written as exp(— iS) where S is the 
eikonal function representing the wave fronts. On a four-dimensional curved 
manifold (M,g), S satisfies the following HJ equation: 

g^dpSdvS = (1) 

g^ v being the inverse metric. To exploit the quaternionic structure of four 
dimensional light cones it is convenient to use a non holonomic base: 

9 = Va b a O b 

Q a = efdx^g iiV = r ]ah efef (2) 
d6 a ^ 

where the greek indices /i run over space-time coordinates, the latin indices 
a = 0, ..,3 run over the tetrad basis and i]^ is the Minkowski tensor on 
tetrad indices. We will denote with n the inverse tetrad. In terms of the 

(a) 

tetrad components eq. (1) reads: 

(4)^) 2 = E(^)^) 2 - ( 3 ) 
i=i 

Equation (3) can be rewritten in the following useful form: 



det O = 0, 9 



s^d^S — \J — le^d^S e^d^S — E^d^S 



(4) 



It is well known that the determinant of a matrix is zero if and only if either 
two rows (two columns) are proportional or one row (one column) vanishes. 

Let us firstly consider the case in which S is a real function. When the first 
row vanishes one obtains: 

e^S + e^S = (5) 
e^S + V^lE^d.S = 0. (6) 

Separating the real and the imaginary part in eq. (6) and rewriting eq. (5) we 
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get 




(7) 
(8) 
(9) 



(4)+ £ (i)Hs = 0. 



The above simple system of equations for S tells us that locally S has a non 
trivial Lie derivative only along the tetrad £/L — so that S can be easily 
found by solving eq. (9) through the method of characteristic. Eqs. (9), (7) and 
(8) could be the starting point of the WKB analysis in gravitational lensing, 
but usually one only considers the characteristics of the above equations, that 
is, usually only the light rays analysis is performed (see, for example, [19]). It 
is clear that S carries more physical informations then single light rays, so, 
as far as gravitational lensing is concerned, this method to treat HJ equation 
could be useful to obtain global properties of the light beams such as the shear, 
the curvature of the wavefronts and so on. 

When the two rows of 9 are proportional (the function of proportionality k 
being an arbitrary non vanishing function on M) and S is real, the equations 
for S can be written as follows: 



The above equations, with a suitable choice of k, can describe interesting phe- 
nomena such as geodesies wrapping round and round in the (2) — (3) direction. 
It could be rather fruitful to exploit the above method together with an in- 
teresting method developed in [12] [13] to deal with the singularities of light 
cones which can be described in a rather detailed way. Indeed, one of the 
most important topic in which the above scheme could be an important tools 
is holography (see, for two detailed reviews, [5] [1]). In the elegant framework 
of covariant entropy bounds [3] [4] which refined pioneering ideas of Bekenstein 
[2] 't Hooft [21] and Susskind [20], an important role to deduce the bounds on 
the entropy of a given two dimensional spacelike surface £ is played by the 
null hypersurface spanned by the non expanding light rays emanating from 
E which, eventually, terminate at caustics. The analysis of such null hyper- 
surfaces is mainly based on geodetics and Raychaudhuri equations. Any null 
hypersurfaces can be seen as a solution of eq. (1) so it is clear that the method 
here proposed could be very important in a holographic perspective: to solve 
eq. (1) is equivalent to solve the null geodetics equation together with Ray- 
chaudhuri equation since the wave fronts S gives informations both on the 
light rays (through its characteristics) and on the shear and on the expansion 
of the light beams (through the induced geometry of the null hypersurfaces 




(11) 



(10) 
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S = const which can be easily studied once S is given). The method here pro- 
posed, which reduces eq. (1) to a much simpler set of linear partial differential 
equations, seems to be an effective tool to explore holography. 

Now, let us consider the case in which S can be a complex function: S = 
V + \J — 1W. It is fair to say that, in this case, the WKB approximation 
breaks down since, with a complex S, it is not true anymore that we have a 
wave with a slowly varying amplitude A and a rapidly varying phase S = V + 
>/— 1W, being the amplitude multiplied by exp (-W) which is rapidly varying. 
The analysis of this case is rather interesting: firstly, it clarifies which is the 
geometrical structure responsible of the breaking of the WKB approximation. 
Moreover, it shed new light on some interesting phenomena occurring in the 
presence of such geometrical structure. A remark is in order here. It is rather 
trivial to construct special complex solutions of eq. (1): given a real solution 
of eq. (1) So, the complex function So + \/— IF (So) is a complex solution of eq. 
(1) for any differentiable real function F(x). Obviously, this kind of complex 
solutions are uninteresting and have no physical meaning. On the contrary, 
the complex solutions we will search for arise only in very special situations 
and are related to geometric structures (which generalize the Cauchy-Riemann 
conditions) signaling interesting physical phenomena. 

When one rows of O vanishes, separating the real and the imaginary part one 
obtains the following linear system of equations for V and W: 

(<o) ± W = 0, (4) ± efo) d,W = (13) 
efadpW = Ts^d.V, efadJV = iefo^V . (14) 

where £/L ±£m is the propagation direction of the geodesic flows and both the 
real and the imaginary part of S have vanishing Lie derivative along the prop- 
agation direction. The non trivial dependence in the transversal directions is 
parametrized by a sort of generalized harmonic functions in the (2) — (3) direc- 
tions since eqs. (14) can be thought as generalized Cauchy-Riemann conditions 
for V and W. To see this, let us consider the case in which 

ef 2) = Cd 2 + Dd 3 , ef 3) = Ed 2 + Fd 3 
CF-DE^ 0. 

With the above tetrads, eqs. (14) read: 

d i W = f 3 i d j V, i,j = 2,3 (15) 
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1 I EF + CD F 2 + C 2 



\ 



CF-DE y_ (£2 + C 2) _ ^ EF + CD j J 



(16) 



thus, we see that an important role in order to have complex solutions is 
played by the tensor defined in eq. (16) being a sort of complex structure 
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in the transversal directions spanned by and e^ 3 y In the following, some 
geometrical conditions on eqs. (14) ensuring the possibility to have non trivial 
complex solutions of HJ equation together with two interesting examples will 
be analyzed. 

In the case in which the two rows are proportional the equations read: 



(1) 



fi)^) 



(17) 
(18) 



kc = A + y/—lB being a non zero complex function on M. Eqs. (17) and (18) 
can be splitted into four real equations as follows: 



L (b) A (a) 
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(19) 



To express the imaginary part of S in terms of the real part, we have to solve 
the linear system (19) for the X^ a y However, det L = since only three rows 
are independent and the first row r (i) is a linear combination of the third 
"^(3) and the fourth ~~r (a) rows ~~r (i) = A~r -\-B~~r In order for the system 
of equations (19) to be compatible, the following relation, which determines 
the real part of S, must hold: 



Y 



(o) 



AY {2) + BY ( 



2B 



(A 2 + B 2 -l) 
(A 2 + B 2 + 1) 



(3) =f 

£ 1i)d,V + 



2A 

(A 2 + B 2 + 1) 



Ay. 



(A 2 + B 2 + l) ( 3 ) M - c (o)-m 



(20) 



By standard linear algebra's consideration, when the above equation is sat- 
isfied, we can discard the first equation of the system (19) and freely assign 
e'Lxd^W, so we will take e^d^W = 0. Eventually, it is convenient, in the re- 
duced system (19), to express the Lie derivative of V along the timelike tetrad 
e^d^V in terms of the Lie derivative along the spacelike ones e^d^V (with 
(i) = 1, 2, 3) through eq. (20). Thus, we arrive at the following system in which 
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only spacelike Lie derivatives are involved: 



= NSWiiW 



iV (i) 



(i) c « L 


2A 
A 2 +B 2 +l 

2B 

\ A 2 +B 2 +l 



r(0 _ 



B 



2A 2 



A B 
B 1 
-A -1 

2AS 



det L« ^ 0, 



J 42 +B 2 + 1 

2AB j _ 



"A 2 +_B 2 +1 



2B 2 
A 2 +B 2 +l J 



The above equation can be formally solved for Xu\ as follows: 



X, 



oiV 



(i) u) M 



(21) 



Now, to have a clearer geometrical picture, let us consider an holonomic co- 
ordinates system in which the spacelike vector fields of the tetrad basis e 



W 

have no timelike components: = 0, V? = 1, .., 3. This, locally, can always be 
achieved. In this coordinates system, the components of the spacelike vector 
fields of the tetrad basis e?U (where we use the first greek letters a, (3, 7 and 
so on, to denote spacelike indices) can be seen as an invertible map between 
spacelike coordinates vector fields d a and the spacelike tetrad vector fields. 
From eq. (21), we can express d a W as follows: 



d n W = T!d«V 



Tl = {e 



(22) 
(23) 



where T% is a tensor field depending only on the metric of the space-time and 
on A and B. Eq. (22) (as well as eq. (15) in the previous case) is related to 
the range of applicability of the standard WKB approximation. 

Let V be a real solution of HJ equation (with suitable choices of k 1 solutions 
of eqs. (10), (11) and (12) are solutions of eq. (20) too), if eq. (22) and eq. (15) 
can be fulfilled with reasonable boundary conditions, then eq. (22) generates 
a new solution with the same real part and an imaginary part representing a 
geometrical magnification or damping factor. Thus, in these cases, the WKB 
approximation cannot be used anymore because the amplitudes of the waves 
travelling in such regions are multiplied by rapidly varying magnification fac- 
tors of the form exp (-W). 

The problem is to determine the geometrical properties of Tp that make this 
possible. Clearly, not any tensor makes eq. (22) (or eq. (15)) consistent: we 



For example k = 1, B = and A = 1. 
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have to impose that W = d 2 a W: 

d 2 a7 W = d 2 a W =► 8, (TPdpV) = d a (Tfd fi V) . (24) 

The solvability of eqs. (24) and the similar compatibility conditions deduced 
from eq. (14) with reasonable boundary conditions puts strong constraint on 
the topology. Let us consider a simple but interesting example, four dimen- 
sional PP-waves (first introduced in [17] [18]): 

g = dt 2 — dz 2 + H(t — z, x, y)d (t — z) 2 — dx 2 — dy 2 , 
(d 2 x + d 2 ) H = 0. 

In this case, eqs. (14), which allows the arising of a complex solution W, 
become: 

d y W = -d x V, d x W = d y V (25) 
This Cauchy-Riemann conditions simply tell us that V is a two-dimensional 
harmonic function and W is the harmonic function conjugated to V. However, 
when one imposes the physical boundary condition that V has to approach to 
a constant in the limit x 2 + y 2 — > oo, one only finds the trivial solution because 
of the maximum principle for harmonic functions. In the case of singular PP- 
waves, (which are naturally coupled to cosmic string along the z— axis and/or 
7— ray bursts propagating in the same direction [9]) the metric and the Ein- 
stein equations reduce to: 

g = dt 2 — dz 2 + H{t — z, x, y)d (t — z) 2 — dx 2 — dy 2 , 

N 



(d 2 x + d 2 ) h = - z)Mt - *)) 



where N is the number of cosmic strings, /Xj is the energy per unit of length 
and (xi(t — z),yi(t — z)) is the position of the i— th cosmic string in the x — y 
plane. Even if eqs. (14) remain formally the same as eqs. (25), now there are 
singularities in the x — y plane corresponding to the positions of the cosmic 
strings. Thus, harmonic functions approaching to a constant in the limit x 2 + 
y 2 — > 00 in this case exist thanks to the singularities in the x — y plane and 
non trivial complex solutions of the HJ equation appear. Another interesting 
example is the Reissner-Nordstrom blackhole: 

t = (1 - ™ + %) # _ U _ ™. + 91) 1 dr2 _ rW . (26) 

In the case in which the 2 and 3 directions are r and (p and 9 = n/2, eqs. (14) 
are: 

TdrW = -djf, d v W = Td r V 




Td 2 + ^d 2 + (rU) r \ :i I / ( I — - ^ ) . (27) 
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If M > \Q\, (in which case the surface gravity is non vanishing) the operator 
appearing in eq. (27) is of mixed type: elliptic outside the two zeros of g u and 
hyperbolic inside. In this case, there is no real continuous piece- wise smooth so- 
lution of eq. (27) approaching to zero for r — > oo (this can be rigorously proved 
using standard techniques in linear partial differential equations). However, if 
M = \Q\ (so that the surface gravity is vanishing), the situation changes and 
real non trivial continuous (but not differentiable on the line r = M) solutions 
appear since the operator in eq. (27) is elliptic degenerate on the line r = M. 
The detailed proof of this fact is trivial but a little bit long. It is possible to 
convince oneself of this by analyzing the Majumdar-Papapetrou solution [14] 
[16] which represents N extreme charged blackholes in equilibrium and is a 
BPS states in supergravity [15]. The metric (26) reads: 

g = U 2 dt 2 - U~ 2 (dx 2 + dy 2 + dz 2 ) 

N 

(d 2 x + d 2 + d 2 ) U = ]TM^, y u Zi ) 

where [x^y^Zi) is the position of the %— th blackhole and Mj its mass equal 
to its charge. In the case in which the 2 and 3 directions are x and y and 
z = Zi, eqs. (14) look like eqs. (25), but there is a singularity in the x — y plane 
corresponding to the positions of the i— th blackhole (xi,yi). Thus, non trivial 
harmonic functions approaching to a constant for x 2 + y 2 — > oo appear. Hence, 
extreme blackholes admit non trivial complex solutions of the HJ equation. 
This characteristic of extreme charged blackholes is indeed related to their 
BPS nature. Non trivial complex solutions appear in the extreme case be- 
cause the two dimensional operator in eq. (27) is elliptic degenerate on the 
line r = M. This happens because the components gu of the metric tensor 
(which approaches to 1 for large r) has only one zero and does not change 
sign and this is a consequence of the fact that extreme charged blackholes ad- 
mit a globally defined timelike Killing vector field 2 . The existence of a global 
timelike Killing vector field is a necessary condition in order to have a BPS 
states in supergravity (see the discussion of Gibbons in [15]). It is a well known 
phenomenon (see, for example, [15] and references therein) that BPS states 
in field theory often carry complex structures similar to the ones in eq. (22) 
and eq. (14). Eventually, when M < \Q\, the operator appearing in eq. (27) 
is elliptic non degenerate and, for any given Dirichlet boundary condition, 
the Cauchy problem admits a unique solution. It is worth to note that the 
problem to understand if a given gravitational field does or does not represent 
an extreme blackhole is reduced to the analysis of Cauchy problems for two 
dimensional linear partial differential equations, this can be a great simplifi- 
cation in more general cases. We learn from this examples that the solvability 
of eqs. (14) and eqs. (22) is able to detect space-time region with not only a 

2 For example, in the Schwarzschild case, the Killing vector dt is timelike outside 
the horizon and spacelike inside. 
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strong gravitational field but also strong electromagnetic fields and/or angular 
momenta. 

The solvability of eqs. (14), (15) and (22) in asymptotically flat space-times 
and, in particular, its character could be related to the surface gravity of the 
horizon in the general case which, unfortunately, is very difficult to analyze. 
The following argument provides an explanation of this relation in the station- 
ary case in which the event horizon is a Killing horizon. It is well known that if 
the surface gravity is non vanishing, the Killing horizon (which can be thought 
as a real solution of HJ equation) is bifurcating while, when it vanishes, there 
is no bifurcation (see, for example, [8]). However, it should be kept in mind 
that, when one says "there is no bifurcation if the surface gravity vanishes" 
one usually means that there is no real bifurcation. In fact, when the surface 
gravity approaches to zero the bifurcation phenomenon does not disappear, 
the new solution simply becomes complex and eq. (24) tells us precisely when 
such complex solutions appear. This explains why the solvability of eqs. (15) 
and (22), through eq. (24), is able to detect Killing horizons with vanishing 
surface gravity, that is, extreme black holes. 

For the above reasons, there is also a relation with the problem of the Cosmic 
Censor conjecture. A Cosmic Censor should not allow gravitational fields 
with generic breaking of the WKB approximation, which characterizes ex- 
treme blackholes, while he should allow only gravitational fields for which eqs. 
(14) and eqs. (22) can not be satisfied, in such a way that gravitational fields 
would not go beyond extremality. These arguments shows that a deeper anal- 
ysis of HJ equation with the above method could shed light on the physical 
constraints needed to make Cosmic Censor conjecture true. 



3 Conclusions 

In this paper, a new method to deal with HJ equation for null hypersurfaces in 
curved four dimensional space-times has been presented. This method, based 
on a convenient form for HJ equation, simplifies the HJ equation without any 
reference to the field equations for the gravitational field, so that it could be 
used even in theories other then general relativity. It allows to find explicitly 
the geometrical conditions on the metric which the gravitational field has to 
satisfies in order to have non trivial complex solutions of the HJ equation. This 
space-times, in which the WKB approximation cannot be used, are of great 
importance since the origin of the breaking of the WKB scheme is related to 
the presence of extreme blackholes. Thus, the analysis of the case in which 
such non trivial complex solutions arise could also shed light on the physical 
conditions to be imposed in order to make the Cosmic Censor conjecture true. 
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